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ABSTRACT 


Elastic  buckling  loads  are  obtained  for  a  family  of  tapered 
columns  using  the  Vianello  iteration  method  as  refined  by  Neumark, 
The  family  of  columns  studied  by  Dinnik  is  generalized  to  include 
columns  not  symmetric  about  mid-length.   An  extensive  table  of 
buckling  factors  is  compiled  with  the  aid  of  a  high  speed  digital 
computer.   Some  of  Dinnik's  results  are  found  to  be  in  error. 
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1.    Introduction 

For  an  initially-straight  slender  column  which  obeys  Hooke's 
law,  the  buckling  load  is  defined  to  be  the  concentric  compressive 
load  under  which  the  straight  configuration  becomes  neutrally 
stable.   The  buckling  load  is  a  function  of  geometry,  material, 
and  end  conditions,  and  it  takes  the  form: 


(1.1) 


where  P  is  the  buckling  load,  m  is  the  buckling  factor,  E  is  the 
modulus  of  elasticity  of  the  column  material,  I  is  the  minimum 
centroidal  moment  of  inertia  of  the  column  cross  sectional  area, 
and  ■#.     is  the  column  length. 

Columns  with  varying  cross-sections  are  often  used  when  the 
prismatic  column  is  inadeguate.   The  varying  cross-section 
column  is  used  in  many  applications  in  order  to  decrease  weight 
or  provide  additional  strength.   Early  studies  of  the  effect  of 
varying  cross-sections  on  buckling  loads  were  done  by  Dinnik  (1) 
and  by  Bairstow  and  Stedman  (2).   Dinnik  presented  a  table  for 
columns  which  tapered  symmetrically  about  a  central  prismatic 
section.   He  based  his  analytical  calculations  on  the  assumption 
of  various  power  laws  for  the  stiffness  in  the  end  portions  of 
the  columns. 

Bairstow  and  Stedman  presented  a  numerical  method  of  solving 
for  buckling  loads  in  columns  with  varying  cross-sections.   This 
method  was  later  improved  by  Ratzersdorf er  (3).   l/ianello  (4) 
presented  a  graphical  method  of  solution  for  buckling  loads,  and 
Morley  (5)  introduced  numerical  integration  into  Viianello's 
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process.   Neumark  (6)  introduced  a  number  of  refinements  to  this 
process  and  summarized  the  entire  procedure.   Approximate 
solutions  based  on  the  energy  method  were  given  by  Bleich  (7), 
and  Boyd  (8)  conducted  both  theoretical  and  experimental  studies 
on  tapered  columns. 

The  present  study  uses  a  numerical  method  of  solution  to 
find  elastic  buckling  loads  for  columns  of  varying  cross-section. 
Specifically,  the  columns  are  tapered,  hinge-ended,  slender, 
initially-straight,  and  concentrically  loaded.   A  large  family 
of  columns  is  studied,  and  Dinnik's  columns,  tapered  symmetrically 
about  mid-length,  are  members  of  this  family.   Dinnik's  data  are 
checked  and  augmented  by  obtaining  additional  points  and  by  pro- 
viding four  significant  figure  results.   The  major  portion  of 
the  results,  houever,  provides  buckling  loads  for  columns  which 
are  not  symmetric  about  midlength.   These  latter  results  are 
believed  to  be  new. 

The  work  can  best  be  described  by  dividing  it  into  three 
phases:   selection  of  a  numerical  method  of  solution;  adaptation 
and  use  of  the  selected  method  in  solving  for  the  buckling  loads; 
and  discussion  and  verification  of  results. 
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2.    Selection  of  a  Numerical  Method 

The  selection  of  a  numerical  method  involved  the  investiga- 
tion of  the  Bairstou-Stedman ,  the  Wianello,  the  Rayleigh  (9), 
and  the  Rayleigh-Ritz  (10)  methods.   Each  of  these  methods  is 
designed  to  find  the  buckling  factor  for  a  varying  cross-section 
column  (11) . 

In  the  discussion  of  the  four  methods,  it  is  assumed  that 
I  ,  the  moment  of  inertia  at  a  longitudinal  distance  x,  is  known 
at  all  points.   It  is  further  assumed  that  the  boundary  conditions 
are  those  of  a  hinged-end  column; 

y  =  D   at   x  -  D  (2.1) 

y  =  D   at   x  =  je  (2.2) 

Both  the  Bairstou-Stedman  and  the  Wianello  methods  utilize 
the  differential  equation  for  the  deflection  curve  of  a  slender 
column  under  a  compressive  load  and  with  small  deflections: 


EL 


<& 


*  c/X2 


Py   -  0 


(2.3) 


uhere  y  is  the  lateral  deflection  at  the  longitudinal  distance 
x.   It  is  convenient  to  modify  this  equation  by  a  form  of 
equation  (1.1) : 


P   - 


-m£h 

1* 


uhere  I„  is  chosen  to  be  the  maximum  moment  of  inertia  in  the 
column.   Equation  (2.3)  then  becomes: 


+  fi(f)y=o 


(2.4) 
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In  the  Bairstou-Stedman  method,  a  value  for  m  is  assumed, 
and  after  the  quantity  ml„/l   is  calculated,  the  fallouing 
approximate  equations  are  used  in  a  numerical  integration: 

/c  =  y?-<  •*•  xyii  (2-5) 

y^i =  y^-i  *  xy/'  (2-7) 

The  primes  indicate  differentiation  uith  respect  to  x,  and  y1*-^-, 

d2 
y"  =  — *>   The  column  is  divided  into  k  increments  of  equal 
dx^- 

length  \,    or i =k X  .   The  subscript  i  indicates  the  longitudinal 

position  and  is  equal  to  kx/j?  ,  i.e.,  when  i  =  D,  x  =  0  and  uhen 

i  =  k,  x  =#.   The  initial  entries  y  and  y"  are  equal  to  zero  in 
'  *o      o 

accordance  uith  the  first  boundary  condition  and  equation  (2.6). 
In  order  to  begin  the  integration  an  initial  value  for  yA    at  x  =  X/2 
is  arbitrarily  assumed.   Then  the  equations  (2.5),  (2.6),  and  (2.7) 
are  successively  employed  until  the  integration  has  proceeded 
along  the  entire  length  of  the  column.   If,  at  the  completion  of 
integration,  the  deflection  corresponds  to  the  second  boundary 
condition,  i.e.,  y.  =  0,  the  correct  value  of  m  uas  assumed.   If 
the  end  deflection  does  not  meet  the  boundary  condition,  neu  values 
of  m  must -be  assumed^  until'  the  second  boundary  condition  is  met. 
In  the  V/ianello  method,  the  initial  deflection  curve  is 

assumed.   Then,  using  equations  (2.6),  (2.7),  and  (2.5),  respec- 

2 
tively,  the  curve  yk  /m  is  obtained  by  numerical  integration. 

By  dividing,  point  by  point,  the  assumed  deflection  curve  by 

2 
this  value  and  multiplying  by  the  knoun  value  k  ,  a  value  of  m 

for  each  point  is  obtained.   If  all  values  of  m  agree,  the  correct 

deflection  curve  was  assumed.   If  they  do  not  agree,  the  curve 
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yk  /m  is  assumed  to  be  the  deflection  curve,  and  the  process  is 
repeated  until  all  values  of  m  are  the  same.   Since  each  assumed 
deflection  curve  is  known  to  be  more  accurate  than  the  preceding 
one  (6),  the  Wianello  process  will  eventually  converge  to  a 
stable  value  of  m. 

The  Rayleigh  method  makes  direct  use  of  the  fact  that  a 
column  is  in  a  state  of  neutral  equilibrium  when  it  is  subjected 
to  the  buckling  load.   This  can  be  expressed  by  the  condition 
T  =  V .   T  is  the  work  done  by  external  loads  as  the  column  is 
displaced  from  a  straight  to  a  bent  configuration  and  \l    is  the 
strain  energy  of  bending  induced  in  the  column  as  it  is  displaced, 
Additionally,  T  =  Pv,  where  v  is  the  decrease  in  distance  between 
the  ends  of  the  column  as  the  column  assumes  the  deflection  curve 


y.   For  small  displacements,  v  is  approximated  by 
i 


2  /  iy')      c/k}  and  so 

T-iPflyfc/K 

0       t 

\i   may  be  expressed  as/ /V  cl*.       where  M  is  the  bending  moment. 

Since  M  =  Py, 

/ 

V    2  r  J   ET* 


Equating  T  and  W  and  solving  for  P, 


or 


■  €7, 


^n  = 


/(/'/ 


'^L   -  & 


a    1 


A  vzc/k      & 


(2. a) 
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The  solution  for  m  consists  of  assuming  a  deflection  curve  y(x) 
which  satisfies  the  boundary  conditions,  equations  (2.1)  and 
(2.2),  and  integrating  to  find  M  and  D.   It  is  generally 
impractical  to  evaluate  D  analytically;  a  numerical  integration 
procedure  is  more  convenient.   An  accurate  value  of  m  can  be 
obtained  by  applying  this  procedure  to  a  number  of  deflection 
curves.   Of  these,  the  correct  or  most  accurate  deflection  curve 
yields  the  minimum  value  of  m. 

The  main  advantages  of  this  system  are  speed  and  simplicity 
of  calculation.   The  choice  of  a  family  of  deflection  curves 
which  can  be  varied  to  produce  the  correct  deflection  curve  is 
difficult,  however,  and  becomes  even  more  so  with  the  intro- 
duction of  the  non-symmetric  column. 

Ritz  provided  an  extension  of  the  Rayleigh  method  which 
eliminates  the  previous  objection.   He  employed  equation  (2.8) 
in  devising  a  method  to  adjust  systematically  the  parameters  in 
an  assumed  deflection  curve  to  yield  a  minimum  value  of  m. 
The  method,  however,  employs  a  comparatively  large  number  of 
separate  numerical  integrations  which  would  become  quite  cumber- 
some in  solutions  for  m  in  non-symmetric  columns. 

The  investigation  of  the  four  methods  was  not  intended  to 
be  exhaustive.   Since  each  method  will  .yield  a  solution,  the 
selection  of  the  method  to  be  used  depended  only  upon  the 
relative  adaptability  of  the  methods  to  the  present  study.   The 
Wianello  method  with  its  automatic  convergence  was  judged  to  be 
most  suitable  to  the  present  study  and  is  particularly  adaptable 
to  non-symmetric  columns.   The  procedure  outlined  by  Mewmark 
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and  the  refinements  which  he  introduced  far  improved  accuracy  were 
used  in  the  solutions.   These  refinements  are  presented  in 
Appendix  A. 
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3.    Solution  for  the  Buckling  Factor 

The  family  of  columns  studied  is  an  extension  of  that  con- 
sidered by  Dinnik.   The  form  is  shoun  in  Fig.  1.   The  moment 
of  inertia  of  crdss-sectional  area  at  distance  x  from  the  left- 
hand  end  is  given  by: 


jl  -n/jt+i//\" 


O  *  x/f  *  _?// 


/* 


h 


-  /.o 


*//*  x/J  ^  (9/jf+*/f) 


A 
I* 


e/4  -*// 


€//f  *-e/t4 


(j/j  +  *//)  *  */4  ^  /  o 


Dinnik  restricted  his  study  to  columns  symmetric  about  midlength, 

i 
i.e.,  b  =  e,  g  =  d,  c  =j?/2.   In  the  present  study,  non-symmetric 

columns  are  also  included.   The  shape  of  an  individual  column  is 

uniquely  determined  by  specifying  I,/I„,  n,  ■a/Jl     ,  and  c//  .   It 

is  easily  shoun  that  the  other  parameters  are  given  by: 

5?//  =  c/£  -    */zf 

c//J?  =    /  -   j/j?   -  <Z/J 


6/1  - 


A 


/-  K 


)3/ji 


/t~ 


K 


/-  K 


</// 


'/» 


where         K  =  (l,/l&)V% 

It  is  of  interest  to  note  that  the  postulated  variation  of 
I  is  attainable  through  linear  taper  for  n  =  1,  2,  3,  and  4. 
Thus:   n  =  1  corresponds  to  a  column  of  rectangular  cross-section 
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CO 
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with  constant  depth  and  linearly  varying  widthj  n  =  2  is  approxi- 
mated by  a  built-up  column  consisting  of  four  angles  connected 
by  diagonals;  n  =  3  corresponds  to  a  rectangular  cross-section 
column  with  constant  width  and  linearly  varying  depth;  and  n  =  k 
is  represented  by  a  solid  truncated  cone  or  pyramid. 

Since  the  liianello  method  automatically  converges,  the 
assumed  shape  of  the  initial  deflection  curve  is  not  critical. 
However,  the  choice  deserves  some  attention  because  the  number 
of  iterations  required  depends  upon  the  accuracy  of  the  starting 
assumption.   The  boundary  conditions  for  a  hinge-ended  column 
require  zero  deflection  at  the  two  ends.   Uhen  1/1      is  unity, 
the  tapered  column  degenerates  into  a  prismatic  column  with  a 
sinusoidal  deflection  curve.   From  preliminary  calculations  and 
observations  of  the  symmetric  columns,  it  was  noted  that  as 
1-,/Ip,  becomes  quite  small,  the  deflection  curve  becomes  approxi- 
mately parabolic.   Uith  these  factors  in  mind,  the  initial 
deflection  curve  is  assumed  to  be: 

y  -  (/-  uz)Ll-I,/lt)  +  l,/l,  s!r{^)  (3.1) 

where  u  =  1  -  2x  /Jt .   Since  the  amount  of  variation  in  the 
deflection  curve  caused  by  the  introduction  of  a  non-symmetric 
column  was  unknown,  the  variation  is  ignored  in  equation  (3.1). 

The  general  solution  used  to  find  the  critical  buckling 
factor  can  be  best  explained  by  an  example,  Table  I.   A  typical 
non-symmetric  column  is  chosen  with  n  =  2,  c/£    =   D.25,  a/^    =   0.0, 
and  1-,/Irp  =  0.5.   For  this  example,  the  length  is  divided  into 
8  segments,  i.e.,  A  =  0.125   and  k  =  8. 
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The  computed  values  of  I  /In  are  entered  on  Line  2.   Lines 

x  2 

3,  k,    and  5  are  self-explanatory.   Line  6  is  the  assumed  value 
of  the  deflection  curve  of  the  column  when  subjected  to  the 
buckling  load.   Line  7  was  obtained  by  employing  equation  (2.6). 


#"=  j&i/TJy* 


(±ht 


(2.6) 


Line  7  is  the  negative  of  the  quotient  of  the  corresponding 
values  on  Line  6  by  those  on  Line  2.   Line  8  is  the  result  of 
altering  Line  7  with  a  refinement  by  IMewmark.   This  refinement 
is  explained  in  Appendix  A.   In  a  symmetric  column,  the  slope 
(y')  is  zero  at  midlength.   For  a  non-symmetric  column,  this 
condition  is  not  met,  but  it  is  nevertheless  used  to  obtain  an 
initial  value  for  Line  9.   Using  this  assumption,  the  initial 
value  on  Line  9  at  x  =  Q.kbS.        is  set  equal  to  one  half  the 
negative  of  the  entry  at  x  =  0.5-£   on  Line  8.   The  entries  to 
the  left  of  x  =  0.^+5^   are  determined  from  right  to  left  by 
successive  subtraction  of  the  intervening  numbers  of  Line  8. 
Line  10  is  determined  by  employing  equation  (2.5). 


yt  -  X--,  *  */r-i 


(2.5) 


The  initial  value  is  obtained  from  the  boundary  condition  that 
y  =  D  at  x  =  0.   The  following  entries  are  the  cumulative  sum 
of  the  entries  of  Line  9.   The  incorrect  choice  of  the  position 
of  zero  slope  and  the  incorrect  assumption  of  a  symmetric 
deflection  curve  lead  to  a  nan  zero  entry  at  x  =  JL      on  Line  ID. 
The  second  boundary  condition  specifies  that  this  value  be  zero. 


2k 


Line  11  provides  a  linear  correction  to  Line  ID  enabling  the 
boundary  condition  to  be  met.   The  entries  of  Line  11  are  com- 
puted by  multiplying  the  negative  of  the  x  =  Jt      entry  of  Line 
ID  by  the  respective  entries  of  Line  1.   The  entries  of  Line 
12  are  the  sum  of  the  corresponding  values  of  Lines  ID  and  11. 
Finally,  the  values  of  m  on  Line  13  are  the  quotient  of  the 
corresponding  entries  on  Line  6  by  those  on  Line  12  multiplied 
by  k2  (64). 

Since  it  is  desired  to  determine  m  to  four  significant 
figures  (three  decimal  places),  the  criterion  for  convergence 
is  that  the  maximum  and  minimum  values  of  m  on  Line  13  be 
within  D.DDD5  of  each  other.   If,  as  in  Table  I,  this  has  not 
occurred,  Line  12  becomes  the  neu  assumed  deflection  curve 
(Line  6)  and  the  process  is  repeated  to  convergence.   At  con- 
vergence, the  value  of  m  is  taken  to  be  the  average  of  the 
maximum  and  minimum  entries  of  Line" 13. 

This  general  method  of  solution  uas  programmed  for  the 
DDC  1604  high  speed  digital  computer.   The  factors  a//   , 
c/j?  ,  I  /I  ,  and  n  were  varied  to  cover  the  usable  range 
of  columns  of  this  type,  and  the  resulting  buckling  factors 
are  listed  in  Table  II.   The  factor  n  varies  from  1  to  4  in 
integer  steps.   For  each  n,  c/j?  varies  from  D.DD  to  D.5D  in 
increments  of  D.D5;  it  should  be  noted  that  values  of  c/Jl 
from  D.5  to  l.D  are  merely  repetitions  of  those  from  D.5  to 
D.D.   For  each  c/J£      ,  1-,/Ir,  varies  from  D.D5  to  l.DD  in  incre- 
ments of  D.D5  and  a/^   varies  from  D.D  to  2c/^   in  increments 
of  D.l.   It  should  be  noted  that  a  column  is  not  defined  when 
a/2  Jl     is  greater  than  c//  . 
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k.         Discussion  Df  Results 

In  a  preliminary  investigation  the  factor  k  was  varied 
from  20  to  160  in  the  program  solution  for  a  column  which 
appears  to  have  one  of  the  most  stringent  reguirements  for 
correct  convergence.   It  was  assumed  that  an  adeguate  value 
of  k  for  this  column  is  adeguate  for  the  remainder  of  the 
columns.   The  column  is  symmetric,  i.e.,  c/J>     =  0.5;  and 
n  -   ^+,  L/L  =  0.10,  and  a/^      -   0.0.   The  following  converged 
values  of  m  were  obtained  for  each  k: 

a)  k  -   2D  m  -  ^.81242 

b)  k  =  kU  m  =  k. 81206 

c)  k  =   80  m  =  if. 81203 

d)  k  =  160  m  =  k. 81203 

It  is  evident  for  the  six  significant  figures  shown,  that  when 
k  =  80,  the  process  converges  to  the  proper  value.   The  reguire- 
ment  is  for  only  four  significant  figures,  however,  and  k  =  ^+0 
safely  fulfills  this  reguirement. 

The  final  results  in  Table  II  were  analyzed  by:   verifying 
that  the  solutions  for  prismatic  columns,  i.e.,  a//  =  1.0  or 
L/L  =  1.0,  agree  with  theory;  verifying  that  the  results  agree 
with  the  few  analytical  solutions  for  which  there  is  a  known 
deflection  curve;  and  comparing  the  results  with  those  buckling 
factors  available  from  Dinnik's  solutions. 

For  a  uniform  slender  column  with  hinge  supported  ends,  the 

2 
buckling  factor  m  is  egual  to  Tf     .   The  buckling  factors  in 

2 
Table  II  agree  with  the  value  of  Tf       to  four  significant  figures 

for  all  columns  where  either  a/ X        or  L/L  is  unity. 
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There  are  four  analytical  solutions  which  are  both  per- 
tinent to  this  study  and  are  readily  obtainable: 

(1)  c//  =D.D,   n=2,   a/i  =0.0,   0*1  /I  *  1.0 

(2)  c/f     =0.5,   n=2,   a#  =0.0,   0*1  /I  -^1.0 

(3)  c#  =0.0,   n=4,   a//  =0.0,   0*1/1-1.0 

(4)  c//  =0.5,   n=<+,   a#  =0.0,   0*1  /I  <  1.0 
These  solutions  were  adapted  from  similar  solutions  by 

Bleich  (12)  and  Timoshenko  and  Gere  (13)  and  are  presented  in 
Appendix  B.   A  comparison  of  the  buckling  factors  obtained 
from  them  with  those  of  identical  columns  from  the  numerical 
solution  in  this  study  is  contained  in  Table  III.   The  agree- 
ment between  the  two  is  excellent.   The  two  non-symmetric 
cases  have  a  maximum  deviation  between  analytical  and  numerical 
results  of  0.00002.   This  is  quite  good  since  the  criterion  for 
convergence  in  the  numerical  program  was  that  the  maximum  and 
minimum  values  be  within  0.0005  of  each  other.   The  maximum 
deviation  of  0.00010  in  the  two  symmetric  cases  is  not  as  small, 
but  is  still  well  within  an  acceptable  limit. 

The  buckling  factors  which  are  applicable  from  Dinnik's 
work  were  obtained  from  Timoshenko  and  Gere  (13)  and  are 
listed  in  Table  IV.   Onder  each  of  these  values  is  the  corres- 
ponding value  from  this  study.   There  is  a  small  deviation  in 
these  entries  in  a  large  number  of  cases;  n  =  1  seems  to  be  the 
prime  offender.   The  maximum  difference  is  0.06.   Since  the 
present  study  agrees  so  closely  with  the  exact  solutions,  and 
since  the  same  method  of  solution  is  employed  at  all  times, 
Dinnik's  results  seem  to  be  in  error  by  the  amount  of  deviations 
from  the  present  study. 
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TABLE  III 
ANALYTICAL  CHECK  UF  THF  CRITICAL 


11/12 
.05 
.  1C 
.  15 
.2C 
.25 

.30 
.35 
.40 
.45 

.rO 

.55 
.60 
.65 
.70 
.75 


85 
9C 
9  5 


N  =  2. 


C/L=  ..0 


NUMERICAL 
RESULTS 
2.30237 
3.59827 
4.21168 
4.73361 
5.19807 

5. 
6. 

6. 
6. 
7. 

62218 
01594 
3  8  581 
73624 
J7045 

7. 
7 

7, 

8 

8 

,39085 
.69931 
.99733 
.2  86  ,9 
.56657 

8 
9 
9 
9 

.83959 
.10586 
.36595 
.62039 

BUCKLING  FACTOR 

A/L=C.O 

ANALYTICAL 
RESULTS 
2.60235 
3.59826 
4.21168 
4.73360 
5.19807 

5.62218 
6.01594 
6.38581 
6.73624 
7. j 7045 

7.39085 
7.59932 
7.99733 
8.28608 
8.56656 

8.83959 
9.1G586 
9.36595 
9.62039 


N  =  4  .  •_ 


I  1/12 

.v5 
.10 

.  15 

'>  ~ 
.  t. ' 

.25 

.30 
.35 
.4C 
.45 
.50 

.55 
.60 
.65 
.7C 
.75 

,8C 

.85 

.90 


C/L=  1.0 

NUMERICAL 
RESULTS 
2.20692 
3.121  J  5 
3.82248 
4.41382 
4.9348  0 

5.40581 
5.83894 
6.24209 
6.62073 
6.97836 

7.31949 
7.64496 
7.95713 
8.25751 
8.54733 

8.82764 
9.-?  99  3  3 
9.36313 
9.61971 


A/L=G.O 


ANALYTICAL 
RESULTS 
2.20691 
3.12104 
3.82248 
4.41380 
4.9348^ 

5.40580 
5.83894 
6.24209 
6.62073 
6.97886 

7.31949 
7.64496 
7.9  5713 
8.2575  • 
8.54733 

a»82764 
9.0  9933 
9.36313 
9.61970 
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11/12 

.05 
.IC 
.15 
.20 
.25 

.30 
.35 
.40 
.45 
.50 

.55 
.60 
.65 
.70 
.75 

.80 
.95 
.90 
.95 

FAtiLE    III     (CONTINUED) 

N-2.0  C/L=).5  A/L=0.0 

NUMERICAL  ANALYTICAL 

RESULTS  RESULTS 

4.64476  4.64466 

5.39883  5,39885  - 

5.93412  5.93414 

6.36478  6.36479 

6.73179  6.73187 

7.05536  7.0554  J 

7.34687  7.34688 

7.61359  7.61359 

7.86047  .7.86045 
8.09101  8.09098 

8.30782  8.30779 

8.51291  S*_51288 
8.70784  8.7G781 

8.99386  8.89383 

9.07199  9.07196 

9.243*8  9.24305 

9.40782  9.40781 

9.56683  9.56682 

9.72069  9.2206.J 


N=4.C  C/L=  .5           A/L=0.0 

NUMERICAL  ANALYTICAL 

11/12  RESULTS  RESULTS 

.05  3.81652  3.81648 

.10  4.812)6  4.8121J 

.15  5.49075  5.49078 

.20  6.02C32  6- .12034 

.25  6.46052  6.46061 

.30  6.84050  6.84055 

.35  7.17659  7.17662 

.40  7.47914  7.47913 

.45  7.755'4  7.75503 

.50  8.00926  8.00924 

.55  8.24542  8.2454J 

.6C  8.46626  8.46623 

.65  8.67393  8.67390 

.7C  8.87014  8.87011 

.75  9.05627  9.05624 

.80  9.23345  9.23342 

.85  9.40263  9.40261 

.9C  9.56461  9.5646o 

.95  9.72016  9.72007 
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TABLE  IV 
COMPARISON  WITH  DINNIK'S  VALUES 


n 

a/e 

h/I2 

0 

0.2 

0.4 

0.6 

0.8 

1 

6.48 

7.5B 

8.63 

9.46 

9.82 

6.484 

7.600 

B. 689 

9.476 

9.819 

2 

5.40 

6.67 

8.08 

9.25 

9.79 

5.399 

6.670 

8.088 

9.251 

9.790 

0.1 

3 

5.01 

6.32 

7.84 

9.14 

9.77 

5.007 

6.303 

7.824 

9.144 

9.776 

4 

4.B1 

6.11 

7.68 

9.08 

9.77 

4.812 

6.113 

7.680 

9.083 

9.767 

1 

7.01 

7.99 

8.91 

9.63 

9.82 

7.009 

7.999 

8.916 

9.555 

9.829 

2 

6.37 

7.49 

8.61 

9.44 

9.81 

6.3S5 

7.4B7 

8.610 

9.444 

9.814 

0.2 

3 

6.14 

7.31 

8.49 

9.39 

9.81 

6.136 

7.296 

8.488 

9.399 

9.808 

4 

6.02 

7.20 

8.42 

9.38 

9.80 

6.020 

7.197 

8.424 

9.374 

9.805 

1 

7.B7 

8.59 

9.19 

9.70 

9.84 

7fB70 

8.614 

9.250 

9.668 

9.843 

2 

7.61 

8.42 

9.15 

9.63 

9. 84 

7.614 

8.430 

9.150 

9.598 

9.839 

0.4 

3 

7.52 

B.  38 

9.12 

9.62 

9.84 

7.524 

8.364 

9.113 

9.621 

9.837 

4 

7.48 

8.33 

9.10 

9.62 

9.84 

7.479 

8.331 

9.095 

9.614 

9.836 

Upper  values  in  each  set  of  numbers  are  from  Dinnik's  work; 
those  from  the  present  work  are  directly  below. 
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TABLE    IV    (CONTINUED) 


n 

a/$ 

Vi2 

0 

0.2 

0.4 

0.6 

0.8 

1 

8.61 

9.12 

9.55 

9.76 

9.85 

a.  ecu* 

9.100 

9.500 

9.751 

9.854 

2 

8.51 

9.04 

9.48 

9.74 

9.85 

8.513 

9.040 

9.469 

9.740 

9.853 

D.6 

3 

8.50 

9.02 

9.46 

9.74 

9.85 

8.482 

9.020 

9.458 

9.7  37 

9.852 

4 

8.47 

9.01 

9.45 

9.74 

9.85 

8.466 

9.009 

9.453 

9.735 

9.852 

1 

9.27 

9.54 

9.69 

9.83 

9.86 

9.262 

9.511 

9.701 

9.816 

9.863 

2 

9.24 

9.50 

9.69 

9.82 

9.86 

9.243 

9.499 

9.695 

9.814 

9.862 

D.a 

3 

9.23 

9.50 

9.69 

9.81 

9.86 

9.237 

9.495 

9.693 

9.813 

9.862 

*♦ 

9.23 

9.49 

9.69 

9.81 

9.86 

9.233 

9.493 

9.692 

9.813 

9.862 
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The  results,  as  mentioned  previously,  cower  the  usable 
range  for  columns  of  the  type  under  study.   There  might  be  a 
necessity,  however,  for  obtaining  intermediate  values  within 
this  range  which  are  not  listed  in  Table  II.   Therefore,  some 
consideration  should  be  given  to  the  method  of  interpolation. 
Since  the  data  obtained  form  curves  of  second  order  or  higher, 
an  error  is  introduced  when  employing  linear  interpolation. 
The  amount  of  error  depends  upon  the  position  of  the  needed 
value  in  Table  II.   The  error  increases  with  increasing  n, 
decreasing  a/J>  ,  and  decreasing  I  /I  . 

A  comparison  of  linear  interpolation  with  interpolation 
using  Bessel's  formula  is  given  in  Appendix  D.  Interpolation 
formulas  and  functions  were  obtained  from  Comrie  (14).  For  a 
column  with  c/^  =0.5,  n  =  4,  I^/Io  =  D.23,  and  a/i  =  D.15,  it 
was  found  that  there  is  a  Q.11%  difference  between  linear  and 
second  degree  interpolation. 


76 


BIBLIOGRAPHY 

1.    A.  I\l.  Dinnik,  "Design  of  Columns  of  Varying  Cross-Section,  " 

Transactions  of  the  American  Society  of  Mechanical  Engineers, 
Vol.  Wt    No.  18  (September  30,  1932),   165-71. 

2-    L.  Bairstou  and  E.  W.  Stedman,  "Critical  Loads  for  Long 

Struts  of  Varying  Sections,"  Engineering,  Vol.  98,  (October 
2,  1914),  403-4. 

3.  0.  Ratzersdorf er,  "Determination  of  the  Buckling  Load  of 
Struts  with  Successive  Approximations,"  Royal  Aeronautical 
Society  Journal,  Vol.  47,  (March,  1943),  103-5. 

4.  L.  Vianello,  "Graphische  Ontersuchung  der  Knickf estigkeit 
gerader  Stabe,"  Zeitschrift  des  Vereins  Deutscher  Ingenieure, 
Vol.  42,  (December  24,  1898),  1436. 

5.  Arthur  Morley,  "Critical  Loads  for  Long  Tapering  Struts," 
Engineering,  Vol.  104  (1917),  295-8. 

6.  N.  M.  Newmark,  "Numerical  Procedure  for  Computing  Deflections, 
Moments,  and  Buckling  Loads,"  Transactions  of  the  American 
Society  of  Civil  Engineers,  Vol.  108,  (1943),  1161-1234. 

7.  Friedrich  Bleich,  Theorie  und  Berechnung  der  eisernen  Bfucken. 
Berlin:   Julius  Springer,  1924.  p.  136. 

8.  J.  E.  Boyd,  "Tapered  Struts:   A  Theoretical  and  Experimental 
Investigation,"  Ohio  State  Oniversity  Engineering  Experi- 
mental Station  Bulletin  No.  25  (1923). 

9.  Lord  Rayleigh,  Theory  of  Sound.   Second  edition.  Vol.  1. 
London:   Macmillan,  1926.  pp  111,  287. 

10.  LI.  Ritz,  Gesammelte  bJerke.  Paris:  Gauthier-Villars ,  1911 
p.  265. 

11.  R.  E.  Newton,  Inelastic  Buckling  of  Columns  of  Varying 
Sections,"  Thesis,  the  University  of  Michigan,  1951.  pp  12-30. 

12.  Friedrich  Bleich,  Buckling  Strength  of  Metal  Structures. 
New  York:   McGraw-Hill,  1952.  pp  186-189. 

13.  Stephen  P.  Timoshenko  and  James  M.  Gere,  Theory  of  Elastic 
Stability.  Second  edition.  New  York:  McGraw-Hill,  1961. 
pp  125-130. 

14.  L.  J.  Comrie,  Chamber's  Shorter  Six-Figure  Mathematical 
Tables.  New  York:  Chemical  Publishing  Co.,  Inc.  1954. 
pp  xxii,  348. 


77 


APPENDIX  A 
NEWMARK  REFINEMENTS 
Newmark  (6)  presents  refinements  uhich  improve  the  accuracy 
of  the  Wianello  method.   He  bases  these  refinements  on  the 
assumption  that  y"  is  a  second  degree  curve,  whereas,  in  the 
basic  method,  it  is  assumed  to  be  linear.   He  assumes  that 
this  curve  takes  the  form  of  a  distributed  loading  curve,  Fig. 
2(a),  and  then  solves  for  the  reactions  at  the  end  of  each 
increment  (A).   He  finds  that: 

nL  "h-(^%  -<-  <$  r    -s)  (A.i) 

fft  =  /?t//?t,-|fr/flr^)       <fl-2) 

Eguation  (A. 2)  is  used  in  the  present  study  when  the  curve 
is  smooth.   When  there  is  a  break  in  the  curve,  eguation  (A.I) 
is  used.   The  breaks  in  the  curve  appear  where  the  tapered 
sections  meet  the  prismatic  section  or  each  other.   The  curve 
is  assumed  to  take  the  shape  of  Fig.  2(b)  and  the  modification 
is : 

/ft  =  2?  (^  %    *"  &r    *~6rt  -j  -s,j 

The  reaction  R  is  the  concentrated  angle  change  (  A   y'). 
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APPENDIX  B 
ANALYTICAL  SOLUTIONS 
(a)  Symmetric  Column,  n  =  2,  a//  =  0.0,  0.0  *  I  /I  <  1.0,  c/i  =0.5 
The  column  is  of  the  form  shown  in  Fig.  3,  and: 

The  differential  equation  of  the  deflection  curve  is: 

or,  substituting  for  I, 


^^   +  ^y   =  0 


(B.l) 


where 


_  2  =   ^/* 
£1 


The  general  solution  of  equation  (B.l)  with  two  arbitrary  constants 

A  and  B  is: 


where 


The  boundary  conditions  for  a  symmetric  hinged  end  column  may  be 
taken  as: 

'z\i 


J- 


-4V  ' 


<7  ^    X  =  O 


$~     /    ,  %J   =  O  *t        X   =  Jf/Z. 

Substituting  these  boundary  conditions  into  equation  (B.2) 

results  in  two  homogeneous  equations  for  the  constants  A  and  B 


As;^Jn(ff] 


+•     B     COS 


f'  <"(i) 


2-»'7.0 


Ah     -h  B/z   -  0 


so 


A 

M 

i  \ 


.LA 


,l 


c 
a 

3 

o 
o 


u  o 
•p  • 
0)  o 

(/>     CO 
CO 
V 

u 

3 
U> 

•r-l 
fa 


N«0 


x 
1 
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Solutions  for  A  and  B  exist  only  if: 

tanU  Jfnf-f)*]   -lh'0 


I, 


which  has  an  infinite  number  of  roots  h.   The  smallest  root  h. 
defines  the  critical  buckling  load.   This  can  be  written  as: 

tan  T      =    4 
*  **  (I,/L) 


where 


r  -  i   h,   In   (1,/L) 


Ulith  a  given  ratio  of  I,/I„,  T      can  be  obtained  and  the  root  h   is 


and 


-m  =  (/+4h,l)[l  -(I/I*)*-] 


(b)    Non-Symmetric    Column,    n    =   2,    a/jp     =0.D,    P.D«  L/L^l.D, 
c/i    -   D.D    or    l.D. 

The  column  is  as  shown  in  Fig.  k    and,  as  for  the  preceding 
column,  the  general  solution  is: 

y   -    /J~~[A  Sfo  (hint))    +    B  Cos  (h  in  J)] 


The  boundary  conditions  are 

Z 
'I 


J  -ffl' 


J  -  *0, 


y-0 


4t     X~   I 


These  conditions  yield  the  two  eguations: 

a  sfaih  in  a//j*7  *  b  cos  [h  Aa/zj^o 


Therefore, 


a  -  0 


J/*/7»   Jfn(l,/lJt)i]  =0 
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a 
a 

3 

1-4 

O  O 

u     • 

1—1 


a  o 
>»   • 


0) 

3 
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or 


A  = 


2  7f 


in  a//j 


and 


(c)  Symmetric  Column,  n  =  4,  a//   =  Q.D,  c/f     =0.5,  0.0*  Ii/I  *  l.Q. 
The  column  is  of  the  form  known  in  Fig.  3  and: 


The  differential  equation  is 


EI   dX 


Py  = 


O 


or 


+■   Py  =  o 


substituting   ^  =  //X.* 

j+\_dt*       t   dt 


Py  = 


=  0 


or 


where 


z4    ^ 

£*  = 

+  S^v  -  a 


(B.3) 


Equation  (B.3)  is  a  form  of  Bessel's  differential  equation  and 
has  the  solution: 

y  =  fi  [a,  t<  (st)  +  Bt  y.4  ^;] 


where 


Bk 


SD 


=  zr-Jjrg      A,  cos  (St)    +  B,s/0r?  (St) 


and  the  general  solution  is 


y-  x 


A,   cor/^r)  +    B.S 


's. 


'/nt-7* 


(B.4) 


The  boundary  conditions  are 


y  =  o 


2^  _ 


Substituting  these  boundary  conditions  into  Equation  ( B . 4 )  and 
solving  for  A   and  B„  yields: 


z«n\ZU)  A 


or 


where 


Then, 


/°  "^-/        2  4  Af^ 


777     -    4  fl  *{/,//*)* 

(d)    Non-Symmetric   Column,    n    =   if,    a//      =D.O,    c/£    -   D.D    or    l.D, 

o.a*  i1/i2^i.d. 

The  column  is  of  the  form  shown,  in  Fig.  k    and  the  general 
solution  is  the  same  as  for  the  preceding  column: 


;-^k«4)  +  s*Mi*) 


(B.if) 
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The  boundary  conditions  are: 

**"  =  6  ,  y  =  o         0^     k  =  o 

x*  =  Jj  y  =  o  *f     x  -  4 

Substituting  these  boundary  conditions  into  Equation  (B.^O  and 
solving  for  A   and  B„  yields: 

uhich  has  an  infinite  number  of  roots  7}7f  .   The  smallest  non- 
trivial  root  7f    defines  the  critical  buckling  load.   Then, 

6J  " 


and 


777  =  If 


zfl\i 
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APPENDIX  C 


PROGRAM  LISTING 
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APPENDIX  D 
INTERPOLATION  METHODS 
A  value  of  m  is  desired  for  c/#     =  0.5,  n  =  4,  1-,/Ir,  =  0.23, 
and  a/S.    =   0.15. 
(a)  Second  degree  interpolation 

Bessel's  formula  is  modified  to: 

m  =  mn  +  pA'  +  P(P"1}  (  A"+A,") 

0/2        4  , 

where  it  is  assumed  that  the  third  differences  are  negligible. 

Il/I2    =   D*15'    a/^     =   D#15 


a// 

m 

A 

A 

0.0 

5.491 

594 

0.1 

6.085 

654 

60 

0.2 

6.739 

693 

39 

0.3 

7.432 

6.085   +    0.5    (.654)    -    .062    (.099)    =   6.406 


h/I2 

=   0.20, 

a//     = 

0.15 

aft 

m 

A' 

A" 

0.0 

6.020 

569 

0.1 

6.589 

608 

39 

0.2 

7.197 

624 

16 

0.3 

7.821 

=   6.589    +   0.5    (.608)    -    .062    (.055) 


6.890 


Ix/I2 

a/jf 

0.0 

0.1 
0.2 
0.3 


0.25,    a/jC     =0.15 

A      A 

536 
560 
562 


m 
6.461 


6.997 
7.557 
8.119 


24 
2 


=  6.997  +  0.5  (.560) 


,062  (.026)  =  7.275 
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Vi2 

=   D.3D, 

a/i    = 

0. 

a/J? 

m 

A 

A 

Q.D 

G.841 

5DD 

Q.l 

7.341 

514 

14 

0.2 

7.855 

505 

-9 

D.3 

8.360 

m   =  7.341   +   0.5    (.514)    -   .062    (.005)    =  7.598 


hn2 

=  0.23, 

a#    = 

0.15 

h/lZ 

m 

i 

A" 

0.15 

6.406 

k&k 

0.20 

6.890 

385 

-99 

0.25 

7.275 

323 

-62 

0.30 

7.598 

m 


=  6.890   +   0.6   (.385)    -.060   (.161)    =  7.131 


(b)  Linear  Interpolation 

Il/I2  =  D,2D'  a/^    =  D*15 
m  =  6.589  +  .304  =  6.893 

I,/I2   =  0.25,    a//    =  0.15 

m   =   6.997    +   0.6    (.560)    =   7.277 

I1/I2   =  0.23,    a/P    =  0.15 

m   =  6.893   +   0.6    (.384)    =  7.123 


(c)  Comparison  between  the  two  methods 

7.131  -  7.123 


difference  = 


7.131 


=  0.11% 
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